Charge qubits and limitations of electrostatic quantum gates 
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We investigate the characteristics of purely electrostatic interactions with external gates in con- 
structing full single qubit manipulations. The quantum bit is naturally encoded in the spatial wave 
function of the electron system. Single-electron-transistor arrays based on quantum dots or insulat- 
ing interfaces typically allow for electrostatic controls where the inter-island tunneling is considered 
constant, e.g. determined by the thickness of an insulating layer. A representative array of 3 x 3 
quantum dots with two mobile electrons is analyzed using a Hubbard Hamiltonian and a capacitance 
matrix formalism. Our study shows that it is easy to realize the first quantum gate for single qubit 
operations, but that a second quantum gate only comes at the cost of compromising the low-energy 
two-level system needed to encode the qubit. We use perturbative arguments and the Feshbach 
formalism to show that the compromising of the two-level system is a rather general feature for elec- 
trostatically interacting qubits and is not just related to the specific details of the system chosen. 
We show further that full implementation requires tunable tunneling or external magnetic fields. 

PACS numbers: 03.67.Lx, 02.70.-c, 85.35.Be, 85.35.Gv 



I. INTRODUCTION 

Quantum computation in its binary concept requires a 
set of two different quantum states, a quantum two-level 
system (qu2LS), that realizes the quantum bit physically 
[y. Some physical systems are intrinsically qu2LSs such 
as the spin 1/2 of a fermion or the polarization of a pho- 
ton. Spin 1/2 systems such as some nuclei of atoms or 
electrons are considered well-suited also because of the 
comparatively long spin dccohcrence times (~ ms to ns) 
lit which derives from the rather weak interaction 
of the spin with its environment. This long coherence 
also has the price that important processes in qubit op- 
eration are rather slow processes. On the other hand, 
decoherence times for charge based quantum systems are 
orders of magnitude shorter (~ ns) due to the compara- 
tively strong Coulomb interaction [a, El- Yet, all quan- 
tum gates including interacting qubits can be expected 
to scale similarly in time and to allow for fast qubit op- 
erations well below the decoherence limits. Furthermore, 
final readout through single-electron transistors (SET) or 
quantum point contacts (QPC) appears to be "straight- 
forward" (i.e., implementable in principle) 00,0. 

In this paper we discuss quantum bits (qubits) encoded 
in the spatial wave function of electrons embedded in 
condensed matter systems, and the emphasis is placed 
on whether it is possible to realize a qubit based solely 
on charge distribution (charge qubit) and capacitive cou- 
pling. For example, having a set of two quantum dots 
(qudots) close enough so that one electron can tunnel 
back and forth, one may envisage a qubit where the elec- 
tron being on one quantum dot (qudot) represents one 
state, and being on another qudot the other state 0. 
Other examples are the cellular automata setups with a 
2x2 array of qudots with excess electrons |9|, [l(J . These 
proposals are fundamentally based on the variability of 
the tunneling t. Yet there are many physical quantum 



systems where the handle on the tunneling is limited or 
non-existent, such as the case of metallic qudot struc- 
tures where the tunneling barrier is determined by the 
thickness of oxide layers in the structure. 

This then suggests the question of whether there is 
a way that electrostatically controlled logical quantum 
gates (qugates) can realize the necessary single qubit op- 
erations and the tunability of interactions between them. 
In this context, it is important to clearly define what one 
means by a quantum bit encoded in a quantum two-level 
system (qu2LS) and what are the requirements for it. 
The following criteria are established for the usefulness 
of a qu2LS Q: 

Al. The qu2LS should include the ground state of the 
system in the working-range of the (tunable) pa- 
rameters; this clearly facilitates the initialization 
process in an experiment and is much more reli- 
able when compared to a qu2LS completely built 
on excited states. 

A2. The qu2LS should be well separated in energy from 
the remaining states in the Hilbert space. This re- 
duces the influence of the remaining Hilbert space, 
whose interference can be insofar interpreted as a 
source of decoherence and transitions to which may 
result in loss of probability in the primary system. 
This lossy channel for the qu2LS is considered in 
more detail below. 

A3. The qu2LS must interact with a set of external 
gates in order to control single qubit states as well 
as the interaction between them without compro- 
mising the two-level system. 

In the case of charge located on a set of well-defined 
quantum dots (a qudot network) , the electrostatic inter- 
action is clearly able to satisfy point (|A3|) where Coulomb 
blockade (or charging) effects introduce a high energy 
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scale (> lmeV) in typical qudots. These structures lo- 
calize the operating electrons and limit the unwanted fast 
decoherence times due to interaction with the surround- 
ing condensed matter environment. Typical decoherence 
sources such as phonons or unstable impurities in the en- 
vironment are regarded frozen out or static, respectively, 
at the low temperatures required for operation of the 
qubit system. Consequently, the primary source of deco- 
herence in this regime is in fact the reservoir of high- lying 
states, and its coherent interference with the qu2LS will 
be consider explicitly in our description. 

The effect of the controlling electrostatic gates on the 
reservoir of high-lying states must be in the adiabatic 
regime. The time dependent manipulation through gate 
action can be estimated in the following manner: the 
gates are considered to act in clearly specified time win- 
dows in a step like behavior: they are turned on and off 
at will. This switching, however, is always carried out 
with a maximum speed which introduces a characteris- 
tic frequency w switch = 2ir/T switc h, where T switch is the 
switching time itself. In order for the influence of the 
reservoir of higher lying states to be negligible, fajJ SW itch 
must be much smaller than the energy difference to the 
closest coupled states in that bath; thus the switching 
must be done smoothly enough (adiabatically) , so as to 
not admix higher states into the lower qubit states. Co- 
herent quantum operation in the qu2LS, however, de- 
mands the switching to be done faster than 1/5, where 
(5 is the splitting of the qu2LS in question, and incorpo- 
rated in criterion (|A2|) . 

With respect to criteria l|A2(l and (|A3|I . an estimate of 
how much of the wave function may be lost for each gate 
operation can be obtained from the Feshbach formalism 
|14| . For the unperturbed low-energy state manifold V 
(in contrast to the remainder of the space Q), an initial 
state fully contained in V will acquire projections in 
Q due to a gate operation V (assumed instantaneous), 
given by 
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where H = Hq + V and Hpq = PHQ is a projection of 
the Hamiltonian, and with the other projections defined 
similarly (see Eq. JIJ below). Space Q is considered to 
be at least an energy Ao separated from space V and the 
change in the matrix elements due to the gate operation 

V is approximated by the splitting <5 induced by that very 

V in the ground state pair in V; with this, the equation 
above can be estimated as 
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The gate operations considered are a sequence of steps 
in the external parameters, which implies that with ev- 
ery one of these steps a small probability fraction is lost 
from the ground qu2LS to the remaining higher lying 
states, and as such it can be considered as an additional 



channel for decoherence, even if the projection is nearly 
reversible in a gate cycle. Moreover, if the ground state 
can be considered sufficiently isolated (5 <C A), the er- 
ror drops quadratically with the ratio 5/ A, so that if 
6 ~ 0.03A, the probability lost per gate operation would 
be smaller than 0.1 %. However, because the gate opera- 
tions will never be performed instantaneously - smooth- 
ing the transitions so that they take longer than Aq~ 
but are faster than - this is clearly to reduce the 
probability loss to the '<2 reservoir' of excited states. 

Thus with proper adiabatic design of the qugates with 
respect to the higher lying 'reservoir', the way to single 
qubit operations is open. However, how exactly these are 
realized still leaves plenty of possibilities, which one can 
imagine being flexible enough. Here the main emphasis 
is placed on capacitively coupled quantum gates and the 
question of whether they allow the necessary single qubit 
operations. Most surprisingly, the answer will turn out 
negative. Despite the great degree of flexibility in geome- 
try of electrostatic gates and system design, we will show 
below that it is not possible to implement fully opera- 
tional qugates without compromising the robustness of 
the qubits. 



II. THE MODEL SYSTEM 

The model network under consideration is a 3 x 3 ar- 
ray of qudots with a single state per site and spin in- 
cluded. For theoretical purposes the structure is taken 
large enough to illustrate the main physics. For an ex- 
periment in this area, however, it is likely more practical 
to choose an array with fewer dots and gates. For the 
analysis, the 3x3 array is flexible and manageable, and 
is used to illustrate our generalized conclusions. Other 
geometries have been explored and yield similar results. 

The 3x3 array of qudots is sketched together with the 
set of external gates in Fig. QJi. As seen from panel (b) , 
only nearest neighbor capacitances are taken into account 
as well as nearest neighbor tunneling between dots. Over- 
all, four parameters enter the model: the capacitance 
from each dot to either one of the gates (C g = 45aF), the 
nearest- neighbor dot-dot capacitance (Cdd = 45aF), the 
dot self capacitance (Cdo — 45aF) and the nearest neigh- 
bor dot-dot tunneling (t ~ 2fieV). For these values, the 
energy cost for double occupancy (the Hubbard U) be- 
comes U = 1 meV, a typically used value in this context. 
Including a dielectric constant of e = 10, these numbers 
represent a typical design where dots of size 200nm are 
separated by approximately also 200nm. These numbers 
can of course change depending on the detailed geometry 
and shapes of the dots. In that sense, the circles in Fig.GJi 
are just a symbolic representation of the dots. However, 
the capacitance values chosen correspond to relatively 
large nano-structures, and the resulting low energy and 
temperature requirements of few tens of mK may thus 
be lifted to some extent by going to smaller structures. 

The Hamiltonian used to describe this system is of the 
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FIG. 1: Model system -3x3 array: (a) schematic layout: 
circles represent quantum dots and the three horizontal and 
three vertical bars represent the gates which are connected 
to the outside world, (b), same as panel (a) but drawn as 
a capacitor and tunnel junction network, where black circles 
represent the qudots. The box symbols in between dots rep- 
resent capacitive tunnel junctions. 



extended Hubbard type 



-^VijhiUj+^Vihi, 



(2) 



where hi = cicjf + c^jCj| in the standard notation. The 
£(i) a refers to the local energy of the state a on the 
i = {1, . . . , n} identical dots and can be used to account 
for the Zeeman splitting of spins in an external magnetic 
field. For most cases, e a is simply set to zero. The tun- 
neling coefficients t?. are considered independent of the 



spin orientation, thus t 



t 



tij, and unequal to 



zero only between nearest neighbors in the qudot net- 
work. The electrostatic energy in the last two terms of 
Eq. (J2J, i.e. the coefficients and Vi, are derived from 
the (total) capacitance matrix of the system (see Add. 1X1 
for more details). 

Throughout this paper, the electronic system of qudots 
is considered to have, for simplicity, a fixed number of two 
electrons (2e) operating in it. Furthermore, since spin 
fhpprocesses occur on a comparatively long time scale 
U y, 13 , they are neglected and thus the overall spin is 
considered constant. The correct statistics for exchange 
of the two electrons is taken into account by the fermionic 
creation (annihilation) operators cf a (c^). Since spin is 
conserved, the basis for the two spin 1/2 particles (2e) is 
conveniently changed to singlet and triplet states which 
are also spin eigenstates. Fermionic statistics constrains 
the spatial wavefunction of the (triplet) singlet states 
under particle exchange to be (anti)symmetric, respec- 
tively. Henceforth, the spin index can be dropped com- 
pletely, and the distinction between singlet and triplet 
states can be incorporated by using (fermionic) bosonic 
operators for (triplet) singlet states, respectively, with 
the constraint of two particle occupation in the system. 




FIG. 2: Potential landscape for array in Fig. (a) Sin- 
gle particle potential on array with no gate voltages applied 
(V g = 0). (b) Change of single particle potential due to one 
specific set of applied gate voltages. 



III. SINGLE QUBIT QUANTUM GATES 

The 3x3 model introduced in the previous section is 
representative for a network with C^ v symmetry. With 
this geometrical symmetry, several states will have a nat- 
ural 2-fold degeneracy. Having set up the total capaci- 
tance matrix for the system, the single particle potential 
landscape for the array is shown in Fig. for the case of 
no gate voltages applied. The absence of nearest neigh- 
bors on the outer boundary results in an electrostatic po- 
tential well such that a single charge stays preferentially 
in the center of the array. Applying a peculiar pattern 
of gate voltages, the 90° symmetry of the potential can 
be broken which effectively alters the potential on the 
middle outer islands only; Fig. depicts an example of 
such a situation. 



A. Numerical Simulations 

Adding a second charge to the 3x3 array leads to a 
competing effect of the potential well structure in Fig. [5^, 
with Coulomb repulsion and the two charges push each 
other halfway outside (see Fig. [SJi+e) resulting in two 
low-lying classically degenerate states with charges ar- 
ranged either vertically J2Ji) or horizontally JU). The 
eigenspectrum for the 3x3 system with two electrons 
and its dependence on the gate voltage pattern shown 
in Fig. 03 is plotted in Fig. The three different spin 
configurations for triplet states have exactly the same 
eigenspectrum. So looking at one specific triplet spin 
configuration, there is still an exact degeneracy in the 
ground state due to the spatial symmetry when no gate 
voltages are applied. The corresponding spatial config- 
urations are shown in Fig. [SJi-l-e. For the singlet states, 
however, a gap opens up. This gap originates from the 
distinct exchange symmetry in the spatial part of the 
wave function when compared to the triplet states. Tak- 
ing as basis states those shown for the triplet manifold in 
Fig. EJi+e, it is clear that for the singlet subspace they 
mix into their symmetric and antisymmetric combina- 
tions (bonding/antibonding states) near the degeneracy 
point where no gate voltages are applied (V s = 0). One 
can consider the tunneling part of the Hamiltonian as a 
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FIG. 3: (a) Energy level spectrum of the 3x3 system and 
dependence of the symmetry breaking pattern of gate voltages 
in Fig. |2j3. Singlet states js,s z ) = |0,0) are shown in red, 
while triplet states |s, s z ) = \l,m = {+1, 0, —1}) are shown in 
dashed black, (b) and (c) Probability distribution over the 3 x 
3 array of the ground pair (qu2LS) for singlet states. Notice 
equal probability for spin up and spin down (|i/>|| 2 = 1^4. | 2 )- 
(d) and (e) Lowest triplet states. Case chosen (s z = +1) 
has only a spin up component (note, however, that spatial 
probability distribution is the same for the sixfold degenerate 
triplet states). 

perturbation (H = Hq + V (t)), so that the extra minus 
sign for particle exchange in the spatial part of the triplet 
wave functions has an interesting effect: the perturbative 
terms that mix the two basis states effectively cancel out, 
so that the degeneracy stays intact. The underlying rea- 
son, as seen below, is that due to the C± v symmetry for 
each path there also exists a mirrored path where the 
particles are exchanged. However, for the singlet states, 
the perturbative terms all come with the same sign and 
add up. Thus the basis states effectively mix, resulting 
in the singlet states shown in Fig. [3t>+c. The gap that 
opens then for the singlet states as a function of gate volt- 
ages, is properly described as an anticrossing in the level 
spectrum and this immediately opens the path towards 
generating Rabi oscillations in the system characterized 
by different charge configurations. 

The gap between singlet states with no voltages ap- 
plied (V g = 0) can be estimated using the Feshbach 
formalism |l4j which is a well-suited approximation for 
a state space well separated from the remainder of the 
Hilbert space and can be thought of as a perturbative ap- 
proach that builds on path histories in the Hilbert space. 
The lowest order term for weak tunneling t and consistent 
with the numerical data is given as 

5-32-^, (3) 

where Ao is the energy gap to the excited manifold (Ao — 
0.03meV in Fig. 3a). This result can be visualized as four 
hops (i 4 ) needed at the cost of at least Aq for each of the 



three intermediate states (Aq 3 ); further, the prefactor 
gives the number of possible low energy paths to go from 
one basis state to the other including particle exchange 
symmetry. 

A numerical simulation of the state evolution that 
makes use of the anticrossing of the singlet state is shown 
in Fig. The state of the qu2LS is represented by a three 
dimensional vector in the Bloch sphere 1] where the (ini- 
tial) eigenstates for V 5 = are taken as the basis for this 
representation. Note that this is a slightly modified def- 
inition of the Bloch vector, insofar as there is a reservoir 
of higher lying states accessible to the system. Panels 
(a+b) show the path of the Bloch vector as three sequen- 
tial qubit operations were performed: the first and the 
last operation are based on the symmetry breaking gate 
voltage pattern shown in Fig. 03, which rotate the Bloch 
vector around the z axis, the horizontal circular path in 
panel (a). For demonstrational purposes, the second op- 
eration is based on tunable tunneling (an effective a x gate 
in the pseudo-spin space of the qubit, see below). This 
second operation rotates the Bloch vector around the x- 
axis and together with the first operation allows one to 
rotate the Bloch vector anywhere in the Bloch sphere as 
required for single qubit operations. Panel (c) shows the 
evolution of the real space probability distribution during 
the second operation. The two basis states in Fig. [SJi+c 
are nicely rotated into each other over a time consistent 
with the nature of the Rabi oscillations. Typical rise- 
times for the voltage gate that do not mix in higher lying 
states are well below the 1 ps range, while if one were 
to tune the tunneling, the minimal rise-times for adia- 
batic switching would require times in the lOOps range, 
in order to limit the admixing of higher lying states. The 
adiabatic regime considered here is clearly seen in the 
Bloch sphere representation of Fig. 0Jd by observing that 
the length of the Bloch vector is reduced to less than one 
in the intermediate gate operation. On return to the ini- 
tial parameters, this amplitude temporarily lost to the 
bath is regained, however. 



B. Rabi Oscillations and 2nd Quantum Gate 

Qubits are conveniently mapped onto the spin 1/2 for- 
malism using Pauli matrices 0, 0|; the system is de- 
scribed by a pseudospin which can be rotated in 31? space 
by applying perturbations which effectively act as mag- 
netic field along different directions (note that there is 
no real magnetic field and that the real spin of the two 
electron system is taken care of by the singlet and triplet 
states) . An arbitrary single qubit operation thus requires 
the realization of two distinct rotations in the 3D pseudo- 
spinor space. This translates to two linearly independent 
combinations of the three Pauli matrices <7{ x . y ,z} required 
to implement the necessary quantum gates. a z is easily 
implemented using the gate voltages and thus applying 
different potentials to different regions in the qudot ar- 
ray as outlined in the previous section. Further, if one 
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FIG. 4: Coherent manipulation of the singlet state under gate 
action - (a) Evolution of the qubit in the Bloch sphere repre- 
sentation after projection onto the basis of the (initial) eigen- 
spectrum at V g = 0. The Bloch sphere is shown in black, 
and the evolution of the Bloch vector in the qu2LS is shown 
in blue, (b) Same as panel (a), but side view, showing the 
slight size reduction due to the adiabatic interaction with the 
higher lying reservoir of states, (c) Coherent Rabi oscillations 
for the sequence A — > B — » A in panel (b) in the direction in- 
dicated with t — WfieV. The probability distribution in real 
space is shown for the 3x3 array over equally spaced time 
intervals in a total time window of 0.56ns which corresponds 
to one period for this tunneling based action. 
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FIG. 5: Numerical exploration of effective (pseudo-) magnetic 
fields from a random sequence of gate voltages (4096 configu- 
rations for every t value), (a) Energy level splitting between 
the lowest two eigenstates (the qu2LS), 5, as well as the level 
splitting between the 2nd and the 3rd eigenstate, A, shown 
with and without gate voltages applied. A well behaved two 
level system exists for t < 0.5 x 10~ 5 eV, while for larger 
t higher lying states cross over, (b) Sampling the gate volt- 
ages randomly, the minimum and maximum pseudo-magnetic 
fields achieved are recorded (H — al + Ba, and thus B has 
units of energy). Since B x ,min and B x ,max are very similar, 
the difference AB X is shown explicitly by the blue dashed line. 
13z,min/max values are clearly discernible. Note that the B x 
is directly related to the gap in the ground state (So in panel 
(a)) 



considers the tunneling t to be approximately constant, 
like in lithographically grown qudot structures with the 
tunneling determined by oxide layers, one may ask if it 
is possible to obtain the second qugate by applying a pe- 
culiar pattern of only capacitively coupled voltage gates. 
The answer turns out negative in the sense that either 
the second gate (o~ x ) is orders of magnitude weaker than 
the first gate (o~ z ) , or it compromises the two-level ground 
state system such that at least one initially well split off 
eigenstate comes within gap distance to the qubit encod- 
ing subspace . 

As a way to illustrate this result in the 3x3 system 
of Fig. [IJ,, all six gate voltages were sampled randomly 
within their parameter space and over a significant range 
of tunneling coefficient t values. The results are shown 
in Fig. [5j From panel (a), the region of the intact two- 
level system involving the ground state is identified as 
the region where 8 (V) <C A (V) and thus t < 5fieV; 
S = Ei — Eq and A = E2 — E\ are the energy differ- 
ences of the lowest three states. It is then clear that the 
condition S (V) > A (V) makes the assumption of an iso- 
lated qu2LS no longer valid. Panel (b) to the right of 
Fig. [S] shows the effective pseudo- magnetic fields defined 
via the most general effective Hamiltonian in the qu2LS, 
namely H = al + B ■ a, where B = B c s stands for the 
effective equivalent of a magnetic field. B z can clearly be 
turned on and off by the gate voltages and ranges from 
zero to the value limited by the applied voltages. Yet, 
B x is overwhelmingly set by the tunneling t and hardly 



responds to different applied voltages. In other words, 
although we have access to a variety of gate voltages and 
diverse ranges, the fixed t-value is the one that essen- 
tially determines the a x gate. As t cannot be varied, it 
negates the qubit control one needs over the entire Bloch 
sphere. The range of B x values obtained by varying local 
voltages is so narrow that it is hardly visible in Fig. [SJd. 
In order to see the difference between the maximum and 
minimum value of B x achieved, their difference lARcj is 
plotted separately: the first dip in this curve is related to 
a change in sign in AB X , while the second kink is already 
well beyond the two-level regime and originates in other 
higher lying states taking over the ground state. 

The numerical results show that despite having access 
to a large set of voltage gates, the second qugate cannot 
be implemented electrostatically under the assumptions 
of constant tunneling and no real external magnetic field 
[l3j . One may argue, that this happens because of the 
peculiar geometry chosen and that there may be other 
geometries which would respond differently. That this is 
not the case will be shown in the following section. 



IV. ELECTROSTATIC INTERACTIONS WITH 
GATES 

From an analytical point of view, some general state- 
ments can be made on the charge states considered here. 
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From our description of the situation above, we can for- 
mulate the following two points: 

Bl. Encoding the qubit in the charge and thus in the 
spatial wave function, demands that the basis of 
the ground state pair be formed from two spatially 
separated wave functions. 

B2. Implementation of a second qugate via electrostatic 
means with the tunneling t kept constant and with 
no real external magnetic field |l3j | results in com- 
promising of the two-level low-energy system, in- 
validating its use. 

We can in fact demonstrate that these two points are 
true in general, as evidenced by the specific geometry 
above. In order to show this, we will use the following 
statements: 

CI. The ground state of any state (single particle, sin- 
glet or triplet) for a real (not complex) Hamiltonian 
must be nodeless, where for the states with more 
than one particle one must consider the restricted 
space i7 = r i < < . . . only within some unique 
sorting scheme, where r*i points to the location of 
particle i (this restriction is necessary since for ex- 
ample on the overall space, the triplet states have 
an intrinsic node due to the particle exchange sym- 
metry) . 

C2. If a matrix element V \ip 2 ) is ^ for a local 
potential, then for 0i ^ 02 at least one of the two 
wave functions must have a node within the space 
fi, and thus must be split off from the ground state 
itself. 

The argument for statement (|C1|) is similar to one 
found in ^5|- The statement follows from the ob- 
servation that any eigenstate ip(ri,r 2 ) with a node 
within Q has a counterpart 0| which has the same 
energy expectation value f n \ip\ ■ H ■ \ip\ = E = J„ ipHip 
and thus by the variational principle, the ground 
state must be always nodeless or, at least, can be 
chosen as such. Statement (IC2I) is shown as fol- 



lows: since J* n ip$ (n, r 2 ) V (ri, r 2 ) V>2 (n, r 2 ) ^ 
with the local potential V (r/, r 2 '; ri, r 2 ) 
V (fx, rb) 5 (fi — ri) 8 (r 2 ' — r 2 ) there must some region 
in space where both %p\ and 02 are 7^ simultaneously. 
Yet, since 0i and tp2 are orthogonal eigenfunctions, in 
order for (V'iIV^) to be = there must be still another 
region in space with both ip\ and ip 2 unequal to zero but 
with a different sign compared to the first region. Thus 
either 0i or ip2 must switch sign from one region to the 
other. 

With this, statement (|B1(1 follows from the observa- 
tion that for some specific set of parameters (within the 
working-range of the qubit) the ground state is degener- 
ate; thus utilizing statement I|C1|) . both of these ground 
states must be nodeless. Yet, they must be also orthog- 
onal to each other, and so similar to statement l|C2|l . 



01 and "02 can be chosen such that ipi is = where 

02 7^ and vice versa. This is what is meant by spa- 
tially separated wave functions. Furthermore, since this 
ground state pair is supposed to be sufficiently decou- 
pled from the remaining states, this situation may only 
change slightly during gate operations. 

For statement <|B2ll . the Feshbach formalism is em- 
ployed once more. In the matrix representation, the 
Hamiltonian of the qubit system with an isolated sub- 
space with index {1,2} is given by 



H 



e x (V fl ) 


■ ■ ■ H w (t) 


... \ 


£2 (V 9 ) 










Hkk' (t) 


H* lk , (t) 


■■■ £k (V 9 ) 






H* k « (*) : 





H PP (V s ) H PQ (i) 
H+ Q (t) %(t,V s ) 



(4) 



with Hpp the projection of the Hamiltonian onto the 2D 
ground state space where S\ — £2 for V g — with V s 
the set of external gate potentials. In this spatial repre- 
sentation, the potential V g enters only in the diagonal 
of the Hamiltonian and, furthermore, the Hamiltonian is 
diagonal when t = 0. So there is no coupling of Hpp to 
the remaining space for t = since Hpq (0) = 0. From 
this structure of the Hamiltonian, the first qugate (cr z ) 
is easily realized by choosing V 5 such that S\ ^ £2', the 
second qugate, however, must be realized through cou- 
pling to the remaining space. For simplicity but without 
restricting the V g is chosen that leaves £1 and £ 2 

constant or just shifts them together uniformly; the effec- 
tive two-level Hamiltonian created by the Feshbach for- 
malism effectively folds the remaining Hilbert space into 
the reduced Hamiltonian Hpp and thus creates a shift in 
the diagonal, as well as generating the off-diagonal ele- 
ments. The latter terms can be straightforwardly related 
to the <t x which gives the splitting in Eq. J3J). Thus by 
comparison, the effective second qugate for singlet states 
is approximated by 



32- 



t 4 



where A(V g ) = A + Ae(V g ) is the gap with applied V g , 
and |Ae/Ao| <C 1, so that the second gate is approxi- 
mated by 



A^ 



1 - 



3A £ (V s ) 



(5) 



Note that this is a maximum estimate since all relevant 
higher lying energies arc supposed to behave collectively. 
Yet, as seen from Eq. JSJ, this second qugate has a much 
weaker dependence on the gate voltages since it must 
be mediated by the coupling t, while the first qugate is 
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sensitive to the gate voltages as Ae(V g ) directly. For 
the ground state two-level system to be sufficiently ideal 
in the sense of decoupled from the rest of the system, 
it must hold that i/Ao <C 1 and also |Ae/Ao| -C 1 for 
the gate operation to not interfere with the higher lying 
states. The consequence is that the initial splitting for 
the singlet states is small and the effect of the second 
qugate only changes this splitting by a fraction which is 
an order of magnitude smaller. In order to get a sig- 
nificant contribution, the second condition |Ae/A | <C 1 
would have to be lifted but that obviously sacrifices the 
two- level system altogether. This proves statement (IB2|I . 

Local electrostatic interaction of voltage gates with a 
qubit system is therefore not sufficient for a full set of 
single qubit rotations. However, revision of the argu- 
ments brought forward clearly leaves two ways out of this 
dilemma: first, the gap (B x ) is controlled by the tunnel- 
ing. Thus tunable tunneling allows for the second qugate 
needed as is well-known [lOj. Second, the Hamiltonian 
was assumed to be real. The argument of a nodeless 
ground state wave function very much relies on that fact 
since for a real wave function a sign change is only pos- 
sible via a transition through zero while in the complex 
case this is no longer required. With this, statement (|B1|> 
becomes irrelevant, and the freedom on {ipi \ V \tp2) ^ is 
greatly increased. Specifically, an external magnetic field 
which makes the Hamiltonian complex, will in fact not 
increase the initially existent (but constant!) a x qugate 
for singlet states, but it can reduce it to zero. Eventually, 
this is again equivalent to an effectively tunable tunneling 

El. 



for the dot-dot capacitances chosen give energy scales for 
the tunneling operations in the constrained two-level sys- 
tem in the few /xeV range. Capacitances of about 45aF 
are related to dimensions of a few hundred nm. There- 
fore, still possible smaller sizes leave room for increased 
energy scales. 
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APPENDIX A: CAPACITANCE MATRIX AND 
QUDOT INTERACTION 

The total capacitance matrix of the system consisting 
of dots and gates is written as 



Ct, 



Cdot — dot 
Cgate-dot Ci 



C dot— gate 
gate— gate 



C2I C22 



(Al) 

where the matrix has been decomposed into convenient 
block notation which separates the dot-dot and dot-gate 
interactions. For a network of capacitors connecting pairs 
of objects (dots or gates), the individual matrix elements 
are given as 



CONCLUSIONS 



fitat _ ) °2,0 T l^j'ji 

13 \-C(i~j) 



dots, gates 



C(i<-*j') for?: = j, 
i 7^ J- 



(A2) 



The use of spatial wave functions to encode quantum 
bits has been analyzed with respect to capacitive elec- 
trostatic interactions. With emphasis on systems such 
as lithographically grown arrays where the tunneling is 
fixed to a great extent by the thickness of the tunnel bar- 
riers (oxide layers) used, it was shown that with constant 
tunneling t and no external magnetic field present, the 
single qubit operations from a system built this way are 
severely limited. The analysis shows that a full set of 
single qubit operations requires either a tunable tunnel- 
ing or an external magnetic field which makes the wave 
function complex and thus introduces further flexibility. 
A uniform external magnet field applied to one qubit, on 
the other hand, needs to be sufficiently localized with re- 
spect to an ensemble of qubits eventually needed. This 
is a very challenging experimental task. The parameters 



where C^q is the necessary self-capacitance of the z-th 
object and C(i <-> j) is the capacitor value between the 
two objects i and j. Here, the capacitance matrix is ap- 
proximated by a network of nearest neighbor capacitors 
where the specific capacitance values enter as parame- 
ters to the model. With the block matrix notation given 
in Eq. (|A1|) . the coefficients Vij and Vi for the Hubbard 
Hamiltonian in Eq. J5J follow as 



e 2 [CfcV V t 



e [V^-CSnCii 1 ], 



(A3) 



The interaction of the dots to the gates is mediated as 
expected by the off-diagonal block C21 of the total ca- 
pacitance matrix and is linear in the applied set of gate 
voltages V g = {V gil ,V g ^, ■ ■ ■}■ 
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